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t^ ■ Abstract 

a^ 

Ground-state properties of spherical even-even nuclei 14 < Z < 28 and 
D , N = 18, 20, 22 are described in the framework of Relativistic Hartree Bo- 

^H ', goliubov (RHB) theory. The model uses the NL3 effective interaction in the 

O ' mean-field Lagrangian, and describes pairing correlations by the pairing part 

of the finite range Gogny interaction DIS. Binding energies, two-proton sep- 
aration energies, and proton rms radii that result from fully self-consistent 
^ , RHB solutions are compared with experimental data. The model predicts 

\^ \ the location of the proton drip-line. The isospin dependence of the effective 

C^ ' spin-orbit potential is discussed, as well as pairing properties that result from 

the finite range interaction in the pp channel. 



PACS numbers: 21.10.D, 21.10.F, 21. 60. J, 27.30, 27.40 



I. RELATIVISTIC HARTREE BOGOLIUBOV THEORY WITH FINITE RANGE 

PAIRING INTERACTION 

X 

C^ ' The structure of proton-rich nuclei presents many interesting phenomena which are very 

important both for nuclear physics and astrophysics. These nuclei are characterized by exotic 
ground-state decay properties such as direct emission of charged particles and /5-decays with 
large Q-values. The properties of most proton-rich nuclei should also play an important role 
in the process of nucleosynthesis by rapid-proton capture. In addition to decay properties 
(particle emission, /3-decay) , of fundamental importance are studies of atomic masses and 
separation energies, and especially the precise location of proton drip-lines. On the other 
hand, nuclear-structure models can be compared in detailed theoretical studies of nuclei with 
a large proton excess. In particular, for proton-rich nuclei in the sd — f shell (10 < Z < 28) 
predictions of the nuclear shell-model can be compared with results of models that are based 
on the mean-field approach. Shell-model calculations of proton-rich nuclei with 37 < A < 48 
have recently been reported in Ref. [|l[, and the structure of proton-drip line around ^^Ni 
has been investigated in the framework of the self-consistent mean- field theory in Ref. . 
In addition to non-relativistic Hartree-Fock and Hartree-Fock-Bogoliubov models, in Ref. 



0] also the relativistic mean-field model has been used, with pairing properties described in 
the BCS approximation. While this approximation is acceptable for nuclei close to the (5- 
stability line, as we move away from the valley of /3-stable nuclei the ground-state properties 
calculated with the BCS scheme become unreliable. For proton-rich nuclei this problem 
might be less critical than for nuclei at the neutron drip-line, but nevertheless we expect a 
much better description of ground-state properties in a framework which provides a unified 
description of mean-field and pairing correlations. 

In the present study we report the first application of Relativistic Hartree-Bogoliubov 
(RHB) theory to the structure of proton-rich nuclei. Models based on quantum hadrody- 
namics have been extensively applied in calculations of nuclear matter and properties of 
finite nuclei throughout the periodic table. In the self-consistent mean-field approximation, 
detailed calculations have been performed for a variety of nuclear structure phenomena [Q . 
For open shell nuclei pairing correlations have been included in the usual BCS approximation 
scheme, but also more consistently in the Hartree-Bogoliubov framework. RHB presents a 
relativistic extension of the Hartree-Fock-Bogoliubov theory. It was derived in Ref. [Q in 
an attempt to develop a unified framework in which relativistic mean-field and pairing cor- 
relations could be described simultaneously. As in ordinary HFB, the ground state of a 
nucleus |$ > is described as vacuum with respect to independent quasi-particle operators, 
which are defined by a unitary Bogoliubov transformation of the single-nucleon creation and 
annihilation operators. The generalized single-nucleon Hamiltonian contains two average 
potentials: the self-consistent mean-field F which encloses all the long range ph correlations, 
and a pairing field A which sums up the pp correlations. The expectation value of the 
nuclear Hamiltonian (non-relativistic or Dirac) < $|if |$ > is a function of the hermitian 
density matrix p, and the antisymmetric pairing tensor n. The variation of the energy func- 
tional with respect to p and n produces the single quasi-particle Hartree-Fock-Bogoliubov 
equations [§] in the non-relativistic framework. In the relativistic extension [Q the Hartree 
approximation is employed for the self-consistent mean field, and the resulting Relativistic 
Hartree-Bogoliubov (RHB) equations read 

ho-m-X A \ fUk{r)\ „ f Uk(r 



A* -hD + m + Xj\Vk{r)J ^' [Vk{r) ) ' ^^^ 

where ho is the single-nucleon Dirac Hamiltonian (|^) and m is the nucleon mass. The 
chemical potential A has to be determined by the particle number subsidiary condition, in 
order that the expectation value of the particle number operator in the ground state equals 
the number of nucleons. The column vectors denote the quasi-particle wave functions, and 
Ek are the quasi-particle energies. The Dirac Hamiltonian 

Kd = -ia-V + (3{m + g^a)+g^uj'^ + gpTspl + e^^^A'^ (2) 

contains the mean-field potentials of the isoscalar scalar a-meson, the isoscalar vector u- 
meson and the isovector vector p-meson. A^ is the electrostatic potential. The RHB equa- 
tions have to be solved self-consistently, with potentials determined in the mean-field ap- 
proximation from solutions of Klein-Gordon equations 

[-A + ml] a{r) = -g^ p,(r) - g^ a\r) - g; a'{r) (3) 



[-l^ + ml]uj\v) = -g^p,{v) (4) 

[-A + m>0(r) = -(7,p3(r) (5) 

-AA°(r) = ep,(r). (6) 

for the sigina meson, omega meson, rho meson and photon field, respectively. The spatial 
components a;, ps and A vanish due to time reversal symmetry. The equation for the sigma 
meson contains the non- linear a self-interaction terms |p . Because of charge conservation, 
only the 3-component of the isovector rho meson contributes. The source terms in equations 
(I) to ® are sums of bilinear products of baryon amplitudes 

P.(r) = E V;^(r)7V,(r), (7) 

P.(r) = E Vli^nir), (8) 

P3(r) = E Vl{v)r,V,{vl (9) 

Pen.(r) = E V;t(r)l^\4(r). (10) 

where the sums run over all positive energy states. The pairing field A in (|l]) is defined 

A,6(r, r') = ^ E ^-"cdlr, r')«:ed(r, r'). (11) 

where a, 6, c, d denote quantum numbers that specify the single- nucleon states. Vahcdk^-, r') 
are matrix elements of a general two-body pairing interaction, and the pairing tensor is 
defined 

«:,,(r,r')= E f/4(r)V^*(r'). (12) 

In the Relativistic Hartree-Bogoliubov theory pairing correlations result from the one- 
meson exchange (a-, lo- and p-mesons) ^. However, if for the pairing part of the interaction 
one uses the coupling constants from standard parameter sets of the relativistic mean-field 
model, the resulting pairing correlations are much too strong. The repulsion produced by the 
exchange of vector mesons at short distances results in a pairing gap at the Fermi surface 
that is by a factor three too large. This is not surprising, since in general the effective 
interactions in the particle-hole and particle-particle channels do not have to be identical. 
In a first-order approximation, the effective interaction contained in the mean-field F is a 
G-matrix, the sum over all ladder diagrams. The effective force in the pp channel, i.e. in 
the pairing potential A, should be the K matrix, the sum of all diagrams irreducible in pp- 
direction. However, very little is known about this matrix in the relativistic framework. And 
although the relativistic theory of pairing presents a very active area of research [|^,P] , only 
phenomenological effective forces have been shown to produce reliable results when applied 
to finite nuclei, especially in exotic regions. In the present work we employ a two-body finite 
range interaction of Gogny type [H in the pp channel of RHB: 



ypp^l^ 2) = J2 e-"'"^-'"^)/'^^)' {Wi + BiP" - HiP^ - M,P"P^), (13) 



with the parameters /ij, Wi, Bi, Hi and Mj [i = 1,2). The pairing interaction is a sum 
of two Gaussians with finite range and properly chosen spin and isospin dependence. The 
Gogny force has been very carefully adjusted to reproduce selected global properties of 
spherical nuclei and of nuclear matter. In the pairing channel its basic advantage is the 
finite range, which automatically guarantees a proper cut-off in momentum space. This 
interaction was employed in the RHB calculations of Ref. |10|. For the DIS parameter 



set of the interaction in the pairing channel, the model was applied in the study of several 
isotope chains of spherical Pb, Sn and Zr nuclei. In Refs. |1TT|- [1^] we have used RHB in 



coordinate space with the DIS Gogny interaction to describe properties of hght nuclear 
systems (C, N, O, F, Ne, Na, Mg) with large neutron excess, as well as ground-states of Ni 
(28 < A^ < 50) and and Sn (50 < iV < 82) isotopes. 

The eigensolutions of Eq. (|l]) form a set of orthogonal and normalized single quasi- 
particle states. The corresponding eigenvalues are the single quasi-particle energies. The 
self-consistent iteration procedure is performed in the basis of quasi-particle states. The 
resulting quasi-particle eigenspectrum is then transformed into the canonical basis of single- 
particle states, in which the RHB ground-state takes the BCS form. The transformation 
determines the energies and occupation probabilities of the canonical states. 

For nuclear systems with spherical symmetry the fields cr(r), ci;°(r), p°(r), and A^{r) 
depend only on the radial coordinate r. The nucleon spinors Uk (Vk) in (|l|) are characterized 
by the angular momentum j, its z-projection m, parity vr and the isospin ts = ±^ for neutron 
and proton. The two Dirac spinors Uk{r) and Vfc(r) are defined 

\^Iu{V)\r)\lji^[ti, v?, s) ) 

g{r) and /(r) are radial amplitudes, Xt is the isospin function, and fijim is the tensor product 
of the orbital and spin functions 



m.,™/2 



The two-component functions 



*„W:^(-^n) and *.«:^(-^:^), (16, 

are solutions of the Dirac-Hartree-Bogoliubov equations 

{hoir) -m- A)$c/(r) + / t/rV'^A(r, r')$y(r') = E^uir) 

Jo 

i-hoir) + m + A)$y(r) + / dr'r'^A{r, r')^u{r') = E^v{r) (17) 

Jo 

The self-consistent solution of the Dirac-Hartree-Bogoliubov integro-differential eigenvalue 
equations and Klein-Gordon equations for the meson fields determines the nuclear ground 



state. In Refs. [0,^,0,0 we have used Finite Element Methods in the coordinate space 



discretization of the coupled system of equations. Coordinate space solutions of the RHB 
equations are essential for a correct description of nuclear structure phenomena that origi- 
nate from large spatial extensions of nucleon densities. These include, for example, neutron 
skins and halos in very neutron-rich nuclei. In less exotic nuclei on the neutron-rich side, 
or for proton-rich nuclei, an expansion in a large oscillator basis should provide sufficiently 
accurate solutions ]TB|JT7[] . In particular, proton-rich nuclei are stabilized by the Coulomb 
barrier which tends to localize the proton density in the nuclear interior and thus prevents 
the formation of objects with extreme spatial extension. In the present work we employ the 
procedure of Refs. [|1^,|1^], and solve the Dirac-Hartree-Bogoliubov equations and the equa- 
tions for the meson fields by expanding the nucleon spinors Uk{r) and Vk{r), and the meson 



fields in a basis of spherical harmonic oscillators for N = 20 oscillator shells |T8[ . However, 
in order to verify that our final conclusions do not depend on the method of solution, for 
nuclei at the proton drip-line we have also performed RHB calculations in coordinate space 



15|. In particular, coordinate space solutions have confirmed our predictions for the location 



of the proton-drip line. 

II. GROUND-STATE PROPERTIES OF PROTON-RICH NUCLEI 

In the present application of the Relativistic Hartree Bogoliubov theory we describe 
the ground-state properties of spherical even-even nuclei 14 < Z < 28 and N = 18, 20, 22. 
While for these neutron numbers the nuclei with 14 < Z < 20 are not really very proton-rich, 
nevertheless they will be useful for a comparison of the model calculations with experimental 
data. We are particularly interested in the predictions of the model for the proton-rich nuclei 
in the lf7/2 region. These nuclei have recently been extensively investigated in experiments 
involving fragmentation of ^^Ni [ll9|-p^. The principal motivation of many experimental 
studies in this region is the possible occurrence of the two-proton ground-state radioactivity. 
In particular, the region around ^^Ni is expected to contain nuclei which are two-proton 
emitters. On the other hand, because of the Coulomb barrier at the proton drip-line, the 
emission of a pair of protons may be strongly delayed for nuclei with small negative two- 
proton separation energies. 

The input for our calculations are the coupling constants and masses for the effective 
mean-field Lagrangian, and the effective interaction in the pairing channel. In the analy- 
sis of light neutron-rich nuclei in Refs. [ pT^P^p!^ , as well as in the study of ground-state 



properties of Ni and Sn isotopes [0, we have used the NL3 parameter set for the effective 

mean-field Lagrangian in the ph channel. The effective interaction NL3 has been derived 

23| by adjusting model calculations to bulk properties of a large number of spherical nu- 



clei. Properties calculated with the NL3 effective interaction are found to be in very good 
agreement with experimental data for nuclei at and away from the line of /3-stability. In 
Ref. [^ it has been shown that constrained RMF calculations with the NL3 effective force 
reproduce the excitation energies of superdeformed minima relative to the ground-state in 
^^^Hg and ^^^Pb. In the same work the NL3 interaction was also used for calculations of 
binding energies and deformation parameters of rare-earth nuclei. In the present study for 
the first time we employ the NL3 effective force on the proton-rich side of the /3-stability 



line. In view of the fact that all the results obtained so far indicate that NL3 is probably 
the best effective RMF interaction, the main purpose of the analysis is to study how well 
the properties predicted by the NL3 force compare with experimental data for proton-rich 
nuclei. However, in order to be more specific in our predictions for the exact location of 
the proton drip-line, we will also use two additional standard RMF effective interactions: 



NLl p5[ and NL-SH p6[. These effective forces have been used in many analyses to calculate 
properties of nuclear matter and of finite nuclei, and generally produce very good results for 
nuclei close to the /3-stability line. In particular, the effective interaction NLl was also used 
in the RHB-|-Gogny calculations of Ref. [|T^] 

In most applications of relativistic mean-field theory pairing correlations have been in- 
cluded in the form of a simple BCS approximation, with a monopole pairing force adjusted to 
the experimental odd-even mass differences [Q. For nuclei far from the valley of /3-stability 
this approach becomes unreliable, especially in the calculation of properties that crucially 
depend on the spatial extensions of nucleon densities. The BCS description of the scattering 
of nucleonic pairs from bound states to the positive energy particle continuum produces an 
unphysical component in the nucleon density with the wrong asymptotic behavior |16|,|T^. 
This effect is more pronounced for neutron-rich nuclei, for which the coupling to the particle 
continuum is particularly important. For proton-rich nuclei the Coulomb barrier confines 
the protons in the interior of the nucleus, and therefore the effect of the coupling to the 
continuum is weaker. However, if pairing correlations are described in the unified frame- 
work of the RHB scheme (or HFB in the non- relativistic approach), the nucleon densities 
display a correct asymptotic behavior. The effective interactions that have been used in the 
pairing channel of RHB are the pairing part of the Gogny force and the density-dependent 
delta force. The finite range interaction provides an automatic cut-off of high momentum 
components, while an artificial energy cut-off has to be included in the calculation with 
zero-range forces. On the other hand, the density-dependent interaction can be adjusted 
to produce surface peaked pairing fields, which can be important for a correct description 
of spatial distribution of densities. A fully self-consistent RHB model in coordinate space, 
with a density dependent interaction of zero-range (delta force), has been used to describe 
the two-neutron halo in ^^Li [B^. In the present study we employ the pairing part of the 



Gogny interaction in the pp channel, with the parameter set DIS 0- 
In Fig. 1^ we display the two-proton separation energies 

S2p{Z, N) = Bp{Z, N) - Bp{Z - 2, N) (18) 

for the even-even nuclei 14 < Z < 28 and N = 18, 20, 22. The values that correspond to the 
self-consistent RHB ground-states are compared with experimental data and extrapolated 
values from Ref. [^]. We notice that the theoretical values reproduce in detail the exper- 
imental separation energies, except for ^^Ca and "^^Ti. In order to understand better this 
result, in Table | we compare the calculated total binding energies for the A^ = 18,20,22 
isotones with empirical values. We find that our model results are in very good agreement 
with experimental data when one of the shells (proton or neutron) is closed, or when va- 
lence protons and neutrons occupy different major shells (i.e. bellow and above A^ and/or 
Z = 20). The absolute differences between the calculated and experimental masses are less 
than 2 MeV. The differences are larger when both proton and neutron valence particles 



(holes) occupy the same major shell, and especially for the N = Z nuclei ^^Ar and ^^Ti. 
This seems to be a clear indication that for these nuclei additional correlations should be 
taken into account. In particular, proton-neutron pairing could have a strong influence on 
the masses. Proton-neutron short-range correlations are not included in our model. 

The results should be also compared with recently reported self-consistent mean-field 
calculations of Ref. 0, and with properties of proton-rich nuclei calculated within the frame- 
work of the nuclear shell model [Q. The calculations of Ref. ^ have been performed for 
several mean-field models (Hartree-Fock, Hartree-Fock-Bogoliubov, and relativistic mean- 
field), and for a number of effective interactions. The results systematically predict the 
two-proton drip-line to lie between ^^Cr and ^^Cr, ^^Fe and ^^Fe, and '^^Ni and ^''Ni. Very 
recent studies of proton drip-line nuclei in this region have been performed in experiments 
based on ^^Ni fragmentation on a beryllium target pT],p2|. In Ref. pl| in particular, ev- 



idence has been reported for particle stability of ^°Ni. In the shell-model calculations of 
Ref. |T]] absolute binding energies were evaluated by computing the Coulomb energy shifts 
between mirror nuclei, and adding this shift to the experimentally determined binding en- 
ergy of the neutron rich isotope. The calculated two-proton separation energies predicted 
a proton drip-line in agreement with experimental data and with the mean- field results 0. 
Compared to the results of the present study, the shell-model total binding energies are in 
somewhat better agreement with experimental data. However, the two models give almost 
identical values for the extracted two-proton separation energies of the drip-line nuclei. The 
self-consistent RHB NL3-I-D1S two-proton separation energies at the drip-line are also very 
close to the values that result from non-relativistic HFB-|-Gogny (DIS) calculation of Ref. 

i- 

By using a fully microscopic and self-consistent model for the calculation of binding 

energies, we have the possibility to analyze in detail the single-proton levels. In Figs. 0, 

^, § we display the proton single-particle energies in the canonical basis as functions of 

proton number for the N = 18, 20, 22 isotones, respectively. The thick solid lines denote 

the corresponding Fermi levels. Although the proton levels do not change much with Z, we 

observe a consistent decrease in the energy splitting between the spin-orbit partners ld5/2 

- ld3/2 and 2p3/2 - 2pi/2 with increasing proton number. We will show that this decrease 

results from the reduction of the spin-orbit term of the effective potential [|14[. The lf7/2 

orbital is unbound for all A^ = 18 isotones, and is very slightly bound for A^ = 20. The 

Fermi level displays a sharp increase with Z for all three isotone chains. In principle, a 

positive value of A should indicate which nuclei are beyond the proton drip-line, i.e. which 

nuclei are ground-state proton emitters. In particular, for ^^Cr, "^^Fe and ^°Ni we find A > 0. 

This is somewhat surprising, since for ^^Fe and ^''Ni the calculated two-proton separation 

energies are positive. We have performed RHB calculations with the effective interactions 

NLl and NLSH, but also for these forces the Fermi level is positive for ^^Cr, ^^Fe and 

^°Ni. For these three nuclei we have also verified the results by performing coordinate space 

RHB calculations. The results are practically identical to those obtained with the oscillator 

expansion method; the Fermi levels for these three nuclei have positive values. Therefore 

it appears that there are cases at the drip-line for which the definition of the two-particle 

separation energy ( p!8[ ) does not correspond to the physical interpretation of the chemical 

potential. 

In Fig. ^ we show the self-consistent ground-state proton densities for the A^ = 20 



isotones. The density profiles display shell effects in the bulk and a gradual increase of 
proton radii. In the insert of Fig. ^ we include the corresponding values for the surface 
thickness and diffuseness parameter. The surface thickness t is defined to be the change in 
radius required to reduce p{r)/ p^ from 0.9 to 0.1 (po is the maximal value of the neutron 
density; because of shell effects we could not use for po the density in the center of the 
nucleus). The diffuseness parameter a is determined by fitting the neutron density profiles 
to the Fermi distribution 

p{r) = Po (l + ea;p(^^— ^)) , (19) 

where Rq is the half-density radius. In going away from the valley of /3-stable nuclei, generally 
the proton surface thickness increases and the surface becomes more diffuse. However, while 
t increases from Si to Ni, the diffuseness parameter a has a maximum aX Z = 20. It appears 
that, as protons fill the lf7/2 orbital, the proton surface becomes slightly less diffuse. This 
could be due to the stronger influence of the Coulomb barrier. In |^ we display the self- 
consistent proton potentials for the A^ = 20 isotones, and in the insert the details of the 
potentials in the region of the Coulomb barrier. We notice how the Coulomb barrier increases 
from 3 MeV for ^'^Si, to 6 MeV in ^^Ni. We include also ^^Ni in our figures for the A^ = 20 
isotones, although this nucleus is not particle stable in our calculations. 

In Fig. ^ we display the proton rms radii for A^ = 18, 20, 22 isotones, respectively. 
The calculated values are compared with experimental data for proton radii from Ref. |^ . 



Except for ^^Si, we find an excellent agreement between experimental data and values cal- 
culated with the NL3 effective force with the DIS Gogny interaction in the pairing channel. 
The model predicts a uniform increase of rms radii with the number of protons. 

In an analysis of ground-state properties of light neutron-rich nuclei |T4], we have shown 



that the relativistic mean-field model predicts a strong isospin dependence of the effective 
spin-orbit potential. With the increase of the number of neutrons the effective spin-orbit 
interaction becomes weaker and the magnitude of the spin-orbit term in the single nucleon 
potential is significantly reduced. This results in a reduction of the energy splittings for spin- 
orbit partners. The reduction in the surface region was found to be as large as ~ 40% for Ne 



isotopes at the drip-line. In Ref. [|1^ similar results were found for the Ni and Sn isotopes. 
The spin-orbit potential originates from the addition of two large fields: the field of the 
vector mesons (short range repulsion), and the scalar field of the sigma meson (intermediate 
attraction). In the first order approximation, and assuming spherical symmetry, the spin 
orbit term can be written as 

Vs.o. = -^VUr), (20) 

r or 

where Vis is the spin-orbit potential [Q 

777 

Vis = (V-S). (21) 

^eff 

V and S denote the repulsive vector and the attractive scalar potentials, respectively, me// 
is the effective mass 



m,ff = m-^iV-S). (22) 

Using the vector and scalar potentials from the NL3 self-consistent ground-state solutions, 
we have computed from ( PO]) - (^) the spin-orbit terms for the N = 20 isotones. They 
are shown in Fig. ^ as function of the radial distance from the center of the nucleus. The 
magnitude of the spin-orbit term Vs.o. decreases as we add more protons, i.e. as we move 
away from /^-stable nuclei. /^From ^^Si to ^^Ni, the reduction is ~ 20% in the surface region. 
The minimum of Vs^. is also shifted outwards, and this reflects the larger spatial extension of 
the proton densities. However, we note that the reduction of Vg.o. for protons is considerably 
smaller than the one calculated for neutrons in Refs. 0,|l3 (~ 35 — 40%). 

The properties of the finite-range and density independent pairing interaction are illus- 
trated in Figs. ^ and |10[ In Fig. |^ we plot the average values of the proton canonical 



pairing gaps Anij as functions of canonical single-particle energies. The gaps are displayed 
for canonical states that correspond to the self-consistent ground state of ^^Cr. Anij are the 
diagonal matrix elements of the pairing part of the RHB single-nucleon Hamiltonian in the 
canonical basis. Although not completely equivalent, Anij corresponds to the pairing gap 
in BCS theory. A very detailed discussion of HFB equations in the canonical basis can be 
found in Ref. fl^. The pairing gaps have relatively large values for deep-hole states. This 



is related to the volume character of the Gogny interaction in the pairing channel. The 
average value at the Fermi surface is between 1.5 and 2 MeV, and Anij slowly decrease for 



canonical states in the single-proton continuum. In Fig. |10| we display the averages of the 
proton pairing gaps for occupied canonical states 

< A, >= ^y^;;f"^ (23) 

where Vnij are the occupation probabilities. The values of < Ap > for the N = 22 isotones 
are plotted as function of the proton number. The average proton gap increases to almost 3 
MeV for ^^Ar, then the pairing correlations disappear at shell closure Z = 20. For the li7/2 
orbital the value of < Ap > is ?s 2.5 MeV. 

In conclusion, this study presents the first application of the Relativistic Hartree Bogoli- 
ubov theory to the description of ground-state properties of proton-rich nuclei. A detailed 
analysis of spherical even-even nuclei with 14 < Z < 28 and A^ = 18, 20, 22 has been per- 
formed. The NL3 parameter set has been used for the effective mean-field Lagrangian in 
the ph channel, and pairing correlations have been described by the finite range Gogny in- 
teraction DIS. In a comparison with available experimental data it has been shown that the 
NL3 + Gogny DIS effective interaction provides a very good description of binding energies, 
two-proton separation energies and proton rms radii. Model predictions for the proton drip- 
line agree with recently reported calculations in the framework of the nuclear shell-model 
and with results of non-relativistic HF and HFB studies. For isotone chains we have also 
discussed the predicted reduction of the effective spin-orbit potential with the increase of 
the number of protons, as well as the resulting energy splittings between spin-orbit partners 
and modifications of surface properties. 
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TABLES 

TABLE L Comparison between calculated and empirical binding energies. All values are in 
units of MeV; empirical values are displayed in parentheses. 



32Si 


269.02 (27L41) 


40 Ar 


343.97 (343.81) 


44Cr 


351.65 (349.99) 


34Si 


284.42 (283.43) 


38Ca 


313.11 (313.04) 


46Cr 


380.19 (381.98) 


36Si 


293.08 (292.02) 


40 Ca 


341.99 (342.05) 


44Fe 


312.07 (-) 


34s 


288.10 (29L84) 


42Ca 


362.95 (361.90) 


46Fe 


352.25 (350.20) 


36s 


307.98 (308.71) 


40Ti 


315.39 (314.49) 


48Fe 


384.42 (385.19) 


38s 


320.77 (321.05) 


42Ti 


348.35 (346.91) 


46 Ni 


306.72 (-) 


36 Ar 


302.52 (306.71) 


44Ti 


373.15 (375.47) 


48Ni 


349.92 (-) 


38 Ar 


327.34 (327.06) 


42 Cr 


314.94 (314.20) 


50Ni 


385.52 (385.50) 
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FIGURES 

FIG. 1. Comparison between RHB/NL3 and experimental two-proton separation energies 
for N = 18, 20, 22 isotones. Black symbols denote empirical values; lines connect symbols which 
correspond to calculated values. 

FIG. 2. The proton single-particle levels for the A^ = 18 isotones. Solid lines denote the neutron 
Fermi level. The energies in the canonical basis correspond to ground-state solutions calculated 
with the NL3 effective force of the mean-field Lagrangian. The parameter set DIS is used for the 
finite range Gogny-type interaction in the pairing channel. 

FIG. 3. Same as in Fig. 2, but for the A^ = 20 isotones. 

FIG. 4. Same as in Fig. 2, but for the N = 22 isotones. 

FIG. 5. Self-consistent RHB single-proton density distributions for the A^ = 20 isotones, 
calculated with the NL3 effective interaction. 

FIG. 6. Self-consistent proton potentials for the A^ = 20 isotones. In the insert the details of 
the Coulomb barriers are shown. 

FIG. 7. Calculated proton rms radii for A^ = 18, 20, 22 isotones compared with experimental 
data. 

FIG. 8. Radial dependence of the spin-orbit term of the potential in self-consistent solutions 
for ground-states of the A^ = 20 isotones. 

FIG. 9. Average values of the proton canonical pairing gaps as functions of canonical sin- 
gle-particle energies for states that correspond to the self-consistent ground state of *^^Cr. The 
NL3 parametrization has been used for the mean-field Lagrangian, and the parameter set DIS for 
the pairing interaction. 

FIG. 10. Average proton pairing gaps < Ap > of the N = 22 isotones. 
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